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Abstract 

Signatures of new physics beyond the Standard Model are often characterized by 
large missing transverse energy {$t) produced in association with multiple jets. The 
dominant Standard Model background to such processes comes from gauge-boson pro- 
duction in association with jets. A standard search strategy involves looking for an 
excess in the m e f / distribution, where m e / / =$t+J2jPj an d Pj denotes the transverse 
momentum of the J-th jet. The region of large m e ff is dominated by jet production 
near threshold, giving rise to large Sudakov logarithms that can change the magnitude 
and shape of the m e // distribution. We present an effective theory framework for the 
resummation of such threshold logarithms. We perform an analysis for exclusive jet 
production using the N-jettiness global event shape, which allows theoretical control to 
also be maintained over large logarithms induced by vetoing additional jets. As a first 
step, we give explicit numerical results with next-to- leading- log (NLL) resummation 
for pp — > 7 + 2 jets in the large m e ff region. 



1 Introduction 



The production of gauge bosons in association with jets is a dangerous background to searches 
for new physics at the LHC. When the gauge boson in question is a Z-boson that decays into 
neutrinos, the resulting missing energy plus multi-jet final state is a dominant background 
to searches for supersymmetry (SUSY) and other models with dark-matter candidates, such 
as Little Higgs with T-parity or universal extra dimensions. Searches for SUSY and these 
other theories in the multi-jet plus missing energy channel typically utilize a shape difference 
between signal and background in the tail of the m e ff distribution, where m e // is defined as 
the scalar sum of the missing transverse energies and jet transverse momenta: 



We do not distinguish between the transverse energy and momentum of a jet. The sum over 
jet transverse momenta may be taken either over only the leading 2, 3 or 4 jets, or may be 
taken over all identified jets in the event; ATLAS [I] follows the first of these approaches, 
while CMS utilizes the second [2]. Considering different jet-multiplicity bins may help probe 
non-standard SUSY theories, such as those with compressed spectra [3] • 

Motivated by its importance, numerous theoretical efforts have been devoted to calculat- 
ing precisely gauge bosons produced in association with jets in the Standard Model. Recent 
advances have allowed W, Z + 3 jets [U El E] and even W, Z + 4 jets |8] to be calculated 
to next-to-leading order (NLO) accuracy in perturbative QCD. These calculations substan- 
tially reduce the residual scale uncertainty of the prediction. However, heavy supersymmetric 
states will populate the tail of the m e ff distributions, where fixed-order perturbation theory 
at NLO may miss important corrections appearing at higher orders. In particular, large log- 
arithms of the schematic form ln(l — m 2 e ^/s), with s the center-of-mass energy squared of 
the hadronic collisions, may be induced by the implicit restriction on soft gluons that comes 
from having nearly all the energy go into the leading few jets and $t- Such effects can be 
exacerbated for gluonic initial states, due to the steeply falling gluon distribution at high 
Bjorken x. They would increase for high m e //, potentially mimicking a SUSY signal. The 
resummation of such soft-gluon logarithms following the approach of Refs. [H [TUJ has been 
previously studied for the related single-inclusive jet production j>t distribution [TT] . and 
numerical results for the Tevatron and RHIC were presented in Ref. [12J. The resummation 
effects were found to be moderate but increasing at high p?, generating the aforementioned 
shape difference. In the 2011 data set, both ATLAS and CMS saw events with m e ff near 
2 TeV, and with the higher- luminosity run of 2012, events even closer to machine threshold 
will be observed. A study of threshold resummation of gauge boson plus multi-jet processes 
at the LHC is therefore warranted. 

We begin a study of large logarithmic corrections to gauge boson plus multi-jet pro- 
duction in this manuscript by considering the next-to-leading logarithmic (NLL) threshold 
resummation for 




(1) 



PP ~ ► 7 + 2 jets 



(2) 
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at the LHC. This process represents a first step toward a study of the W/Z + n jet process, 
but is also interesting on its own as a possible calibration process for missing energy plus 
jet backgrounds [EJ. We utilize the Soft-Collinear Effective Theory (SCET) [H HSl [T6l rj7] 
to study the effect of threshold logarithms. A formulation of threshold resummation within 
SCET applicable to multi-jet processes was given in Ref. [18J, and we build on the approach 
outlined there. In addition we incorporate N-jettiness |19j . an event-shape based formalism 
for exclusive jet production, in our analysis. The N-jettiness formalism allows one to veto 
additional jets and maintain theoretical control of the induced logarithmic corrections to all 
orders in perturbation theory. All jet-algorithm dependence is power suppressed. Such an 
event-shape based analysis simplifies higher-order calculations, since these do not depend 
on the jet algorithm up to power corrections. The production of electroweak gauge bosons 
at high transverse momentum, which is inclusive in the recoiling hadronic radiation and 
consequently does not a require jet definition, has been previously studied [20] using SCET. 

We derive a factorization formalism for the production of a gauge boson in association 
with N-jets in the threshold limit that takes the schematic form 

da~H®Ji®---®J N ®S®f®f. (3) 

Here H is the hard function that encodes the physics of the production at the hard scale and 
can be obtained from known fixed-order calculations, Jj denote jet functions which encode the 
effects of collinear radiation within the iV-jets in the final state, 5* describes soft radiation 
inside and outside the N-jets throughout the event, and / denotes the standard parton 
distribution function. The cross-section in Eq. ^ is inclusive in $t and the transverse jet 
momenta subject to the constraint of Eq. Q. The phase space in the region of large m e ff 
is then dominated by configurations corresponding to N hard jets with only soft radiation 
occupying the region between the beam and jet directions. We discuss this point later in the 
text. In the threshold region, factorization is given in terms of the initial state PDFs and 
does not involve beam functions [2T] that can arise in other processes away from threshold 
where energetic collinear radiation is allowed in the beam directions. 

Another potentially large logarithmic correction can arise when the gauge boson is either 
soft or collinear to either the initial state or one of the final-state jets. However, experimental 
analyses often demand that the ratio ^T/ m eff be greater than some minimum value, and 
impose a minimum angular separation between the missing-energy vector and each jet [3]. 
These constraints reduce the effect of such terms, and we consequently do not consider their 
effect here. 

We present numerical results for the 7 TeV LHC production of 7 + 2 jets, and in particular 
study the enhancement of the NLL cross section over the leading-order (LO) result. The 
ratio of the NLL result over the LO one is approximately 1.5 at m e ff ~ 2 TeV, and increases 
further closer to machine threshold. This shape difference demonstrates the importance of 
having threshold logarithms under theoretical control in searches at the LHC. 

Our paper is organized as follows. We review the kinematics relevant to the study 
of 7 + 2 jet production in Section [2j We show here that the phase-space region with a 
large separation between jet transverse momenta gives a suppressed contribution to the 
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high m e ff distribution. In Section [3] the operator basis we use to match QCD onto SCET 
is introduced. We discuss the factorization theorem for the threshold region in Section |4j 
Numerical results are presented in Section [5j Our formalism confers theoretical control over 
logarithms associated with vetoing extra jets in addition to threshold logarithms. We discuss 
potential applications of this feature, as well as open questions and other future directions, 
in Section [6} A detailed description of the hard-matching coefficients and the solutions to the 
renormalization-group evolutions equations for the hard, soft, and jet functions, are given in 
the Appendix. 



2 Kinematics 

In this section we describe the kinematics and relevant degrees of freedom for the production 
of N-jets together with a gauge boson near threshold. We also elaborate on the definition of 
the threshold region which dominates the phase space in the large m e ff region. We closely 
follow the work of Ref. [18] and adapt their formalism to an event-shape analysis using 
N-jettiness [19] instead of traditional jet algorithms. 

Each jet is characterized by its transverse momentum and pseudorapidity rjj. We 
demand a minimum partonic center of mass energy s min defined as 

N 

Smin = (q + ^Pj)*, (4) 

i=l 

where the massless momenta p l j are defined as 

P.J = (Pj cosh r]j, pj, p T j sinh r]j) (5) 

and q is the gauge-boson momentum. The ratios between s m i n and the partonic and hadronic 
center-of-mass energies squared s and s are characterized by the variables z and r: 

z = — — , r = , r < z < 1. (bj 



The limit of hadronic threshold r — > 1 automatically forces the partonic threshold z — > 1. 
However, other dynamical effects such as the steepness of parton luminosities can force the 
partonic threshold limit z — > 1 even away from hadronic threshold [221 [23]- 

In the limit of partonic threshold, the initial-state collinear partons can emit only soft 
gluons. Additional collinear radiation arises from emissions off the final-state hard partons 
that form the jets. Thus, the relevant degrees of freedom in the limit of partonic threshold 
correspond to soft and collinear modes with momentum scalings 

collinear: p c ~ V^A 2 , 1, A), soft: k s ~ V^A 2 , A 2 , A 2 ), (7) 

where A ~ \J\ — z. The momenta are decomposed in terms of light-cone coordinates as 
P = (p + iP~ ,P±)- The collinear modes along each jet are decomposed in terms of light- 
cone coordinates with the spatial components of the light-cone vector aligned along the jet 
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direction. Momentum conservation at the partonic level is given by 

N 

Pl = q + ks + J2Pi, (8) 

i 

where pi denotes the total initial state partonic momentum so that p\ = s. The total final- 
state partonic soft momentum is denoted by k s , while p l c denotes the momentum of the hard 
parton that eventually forms the i-th final-state jet. 

To define N-jettiness [19], the final-state particles are grouped into regions that are as- 
sociated with either one of the jets or the beam directions using a well-defined distance 
measure. The total momentum of a jet then corresponds to the sum of the momenta of 
all the particles grouped into the corresponding jet region. The total N-jettiness for the 
event is defined as 

N+2 

7* = E^ = E 2& ' Pi; p " = Y,Pk e & ■ ^ - & ■ Pk). (9) 

i i k 

The qi denote massless reference vectors along the N jet and two beam directions, Pj is 
the total momentum of all particles grouped into region z, and t 1 n is the contribution of 
region i to the total jettiness r^. The total momentum of the z'-th jet is defined to be Pj, 
corresponding to the total momentum of all particles in the z'-th region. One can also derive 
cross-sections that are differential in the t % n of the various jet and beam regions. 

All the collinear particles associated with the z'-th jet direction are naturally grouped into 
the z'-th region. The total final-state soft momentum can be decomposed as 

N 

k s = Kut + ^2 ki > ( 10 ) 

i 

where fc, is the momentum contribution of the soft radiation to the z'-th jet and k out is the 
total soft momentum that is not grouped with any of the jets. In other words, the momentum 
k out is defined as the total momentum of soft particles that are grouped with one of the two 
beam directions. 

It is useful to decompose the total final-state momentum into two parts. The first part 
corresponds to the minimum momentum needed to create the color-neutral sector with mo- 
mentum and N-jets with massless momenta pj as in Eq. ([5]). This minimum total mo- 
mentum leads to a partonic center-of-mass energy squared given by s min of Eq. Q. The 
remaining part of the total momentum leads to the actual partonic center-of-mass energy 
squared s > s m in- This decomposition can be made explicit by noting that any four momen- 
tum p^ can be written in terms of the massless momentum pj of Eq. ^ as 

jf = p ^ + pV, (li) 

where = (1, 0) andp + = p° — |p|. This follows fromp = p°j + (p°— p°j) andpj = |pj| = |p|. 
Thus, the partonic momentum conservation can be written as 

N N 

P? = 5" + Km + Ylp"i + v " i 5>* c+ + k t)] . (12) 

i i 
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which is equivalent to the condition 



Or N n 

1 - z = PI ■ k out + p? J>f + + (9(A 4 ). (13) 
s L *— ' J 

Thus, the limit of partonic threshold constrains the energy component of soft radiation 
outside the jets (k out ) and the null components of the total jet momenta (p^ + + kf). We note 
that the threshold variable z must be defined carefully in order to avoid the appearance of 
non-global logarithms [24]; the definition above is expected to not contain such terms [18J. 

We note that there can be configurations consistent with Eq. ([I]) where one or more of the 
N-jets becomes soft, so that there is a hierarchy between the different transverse momenta 
Pj. However, the inclusive nature of the observable we consider forces such configurations to 
occupy a small corner of phase space. They are consequently phase-space suppressed. This is 
the same phase-space suppression that makes fully inclusive Drell-Yan processes insensitive 
to special exclusive jet configurations. In our case, we are interested in the tail of the m e ff 
distribution so that m e ff > 1 TeV. Following the experimental studies, we further restrict 
P T j > Pmin with p^in ~ 100 GeV, and impose standard cuts demanding that the photon and 
jets are well separated. Two types of regions can be identified. In the first region, all the jet 
momenta are of the same order so that |pj. | ~ \pj \ . A second type of region can arise where 
one or more of the jets becomes soft so that hierarchies such as |pjj 3> |pjj can arise. This 
second type of region is suppressed for several reasons. First, since p^ in <C Tn e ff, events with 
widely disparate jet momenta populate only a small corner of the full phase space. This can 
be understood by noting that if one of the jets is soft, then the integration measure over the 
corresponding jet transverse momentum scales like (Pmin/ m eff) 2 - Similarly, the integration 
measure over its rapidity is also suppressed in order to maintain the requirement that the 
jet is soft. The allowed phase space of the remaining hard jets is also further restricted to 
the high transverse-momentum region by the requirement that the event have large m e ff, 
which must now be accomplished with a fewer number of states. In other words, since we 
are inclusively integrating over each jet transverse momentum up to the maximum allowed 
valued determined by m e ff, we are not sensitive to special configurations that occupy a 
small corner of phase space. Another source of suppression comes from the hard production 
amplitude. One or more low-pr jets require that the remaining jets and gauge boson have 
very high p T in order to produce a large m e ff. This configuration is suppressed by the steep 
fall-off of the matrix elements at high transverse momentum of the leading jets. 

For illustration we show numerically the suppression of these configurations for the pro- 
cess pp — > 7 + 2 jets in Fig. [TJ In the left panel the average ratio of the sub-leading jet 
Pt over the leading jet px is plotted. This ratio is approximately 0.7 in the high m e // 
region, indicating that the jets have roughly the same pr on average. Further evidence 
is shown in the right panel, where a soft-jet region of phase space is defined by requiring 
100 GeV < p\ < 0.2 x m e ff/2. The contribution of one soft jet and one hard jet to the 
high-m e // distribution is about 5% of the contribution from two hard jets, indicating that 
it contributes only a small amount to the cross section. These numerical results are in fact 
a conservative estimate of the suppression of the one soft-jet region, since a portion of the 
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region with a low-p^ second jet includes a second jet that is energetic but at higher ra- 
pidity. The region of large m e // is therefore dominated by the production a gauge boson 
and N-hard jets near threshold. The dynamics of this configuration is best described by 
N-collinear sectors and soft radiation with threshold kinematics as described in this section, 
and is well-suited to the application of the SCET formalism used here. 




Figure 1: The left panel shows the average ratio of the second jet pr over the pr of the 
leading jet for the ug partonic channel (magenta dotted line), the uu partonic channel (blue 
dashed line), and for the total result (red solid line). The right panel shows the ratio of the 
cross section for one hard jet and one soft jet over the two hard-jet cross section, where the 
soft jet is defined by 100 GeV <p\< 0.2 x m e fj/2. The double-hard region is then defined 
as the difference of the cross section with both transverse momenta greater than 100 GeV, 
minus the cross section with either jet having a pt between 100 GeV and 0.2 x m e ff/2. 
These results use the LO approximation of the cross section. 



3 SCET operator basis 

The process pp — > 7 + 2 jets receives contributions from various partonic channels, which we 
enumerate here. The channels with four participating quarks or antiquarks are given by 

qq' qq'ry, qq' qq'-f, qq' qcfj, (14) 

and the channels with two quarks or antiquarks and two gluons are given by 

qg -> <igi, qg -> qgi, qq -> ggi, gg -> qqi- (15) 

Contributions where the two jets are produced by the decay of an electroweak gauge boson 
are expected to be small, and are neglected. The processes above are mediated by operators 
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in SCET with Wilson coefficients determined by a matching calculation from QCD to SCET. 
The SCET operators can be decomposed in terms of their color structures. There are two 
color structures for the four-quark operators, 



a/37<5 



e. 



m t a xf) (xit a xt), 



(16) 



and three color structures for the operators with two quark fields and two gluon fields: 

©it = {x a 2t a H^xi)A^Al 



1«3 



(X2 t a3 t ai xl) ^M" 3 



(17) 



The Dirac indices are denoted by {a, (3, 7, 5}, with the remaining Greek indices denoting 
Lorentz indices. The amplitudes for the four-quark processes take the form 



1=1 



where e M is the photon polarization vector, 1 g is the matching coefficient that also in- 
cludes the spin structure, and (9^ lS ) denotes the matrix element of the corresponding SCET 
operator. Similarly, the amplitude for the channels with two gluons is given by 

3 



Me 



(19) 



1=1 



The Wilson coefficients in Eqs. (18) and (19) depend on the specific partonic channel in 



Eqs. (14) and (15) under consideration. 



4 Factorization 

In the threshold limit, the process pp — > 7+N jets is characterized by N-narrow jets and only 
soft radiation outside of these jets. The dynamics of such a process can be described in terms 
of collinear degrees of freedom along the jet directions, soft degrees of freedom throughout 
the event, and the initial-state PDFs. The factorization formula takes the schematic form 

da rsj H u ®S JI ®J 1 ®---®J N ®f®f, (20) 

where (g) denotes a convolution structure, Hjj denotes the hard function and Sji denotes 



the soft function with the indices I, J running over the color structure basis in Eqs. (16) 



and (17). The factors of Jj correspond to the collinear jet functions and / denotes the PDF. 



The hard function encodes the physics of the hard partonic interaction, the jet functions 
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describe the dynamics of the collinear momenta in the jet regions, and the soft function 
describes the soft radiation in the event which is either grouped in one of the jet regions or 



outside of the jet regions as in Eq. (10). The momentum scalings of the soft and collinear 
degrees of freedom in the threshold limit are given in Eq. ([7]). The relevant scales in the 
problem can be characterized by 

v 7 ! > v^Vl - z > \Ts{\ - z) > Aqcd, (21) 

so that the hard function, the jet functions, and the soft function are evaluated at the typical 
scales fx h ~ V§, £tj ~ V§y/1 — z, and Lis ~ V§(1 — z), respectively. The optimal choice 
of scales for the most stable resummation can be affected by the shape of the luminosity 
function (22J [23]. This is discussed in the next section. 

In the N-jettiness formalism, another scale that appears is related to or equivalently 
the t 1 n in Eq. The t 1 n correspond to the contribution to the total N-jettiness from region- 
i in the event. For the process pp — > 7+2 jets, this is divided into four distinct regions as 
shown in Fig. [2] Regions 1 and 3 correspond to the two beam directions, while regions 2 
and 4 correspond to the two jet regions. We give a factorization formula that is differential 
in the quantity 

f 2 = r/ 2 +r 2 j4 . (22) 

where r/ 2 and r/ 4 correspond to the contribution to the 2-jettiness from the regions of 
jet-2 and jet-4 respectively. Note that f 2 is distinct from r 2 since it does not include the 
contributions from the soft radiation in regions 1 and 3 which lie outside the two jet regions. 
In the threshold limit, f 2 ~ V^X 2 , corresponding to the size of the jet invariant masses. 
Thus, f 2 has a size that corresponds to the typical soft or residual momenta in the event. 
The factorization formula takes the general form 



da (2tt) 



1 rhln 



2 - 



x5' 2 (qi+ qi+ q r )5 [y-tanh 



dz I dy I ds j2 ds j4 I dk° out / dr 2 j2 s dr 2 J t 



'~j 2 +qj 4 + q )5(2/-tanh ( ^o +q o +q o 



V Vs Vs J V QJ2 <?J 4 / 

xf(x a} fjLs) f(x b ,Li s ) Jj 2 (sj 2 ,li s ;lij) Jj 4 (sj 4 ,li s ;lij) , (23) 



where we have defined k out = k out -pi/\pi\- The cross-section above is fully differential in the 
massless jet momenta qj 2 and qj 4 defined in Eq. ([5]), the photon momentum q, and f 2 . We 
note that c?$j denotes the phase space of the massless particle i, 

= <PP* . (24) 
1 2(2ixfE i v ; 
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Figure 2: The four jettiness regions into which each event is divided, illustrated using the 
QQ ~~ 991 process. The dashed lines indicate the separation between the regions. Regions 1 
and 3 correspond to the beam directions, while regions 2 and 4 correspond to the two jets. 

The initial partonic momentum fractions x a> b are given by 

x a = ^~e v , x b = J^e~ v - (25) 

The arguments q~i of the hard function denote the projection of the total final-state mo- 
mentum onto the direction rij, with % denoting any of the four regions defined above. This 
factorization formula can be easily generalized for processes with more than two jets. As 



shown in the schematic form given in Eq. (20 ), for the case of N-jets the above formula will be 
modified to have an appropriately generalized soft function convoluted with N jet functions. 
The m e ff distribution can be obtained from the factorization formula above by performing 
the integrations over the remaining phase space factors and over T2 after inserting the delta 
function constraint 6(m e ff — qj 2 — qj 4 — q T ). 



The specific form of the hard, jet, and soft functions in Eq. (23) will depend on the 
partonic channel that mediates the process pp — > 7 + 2 jets. For example, for the partonic 
channels where the final-state jets are initiated by quarks, the factorization formula will be 
in terms of quark jet functions. Similarly, processes with jets initiated by gluons will involve 
gluon jet functions. The color representations of the Wilson lines in the soft function will 



also depend on the partonic channel. Similarly, the PDFs in Eq. (23) will correspond to 



that of the initial-state partons in the partonic process. Finally, the hard function must 



9 



be computed separately for each partonic channel, and is given by the spin-summed and 
color-ordered partial amplitude squared. 

The field theoretic definitions of the different types jet and soft functions are listed in the 



Appendix for completeness. For the partonic channels in Eq. (14) and Eq. (15), mediated 



by the operators defined in Eqs. (16) and (TTTJ) respectively, the hard function takes the form 



u e — r, n 1 ^ n* J irf (A \aa' 
n IJ — ~9w' { ^a/3"/5 { ^a'0'Y6'yn2) 



(26) 



and 



U IJ — ~9pp'9p.p,'9vv l{ ^ a R ^a'R' 



(27) 



The Wilson coefficients C^g s and C 1 ^ 11 are defined in Eqs. (18) and (19) and depend on 



the partonic channel being considered. The momenta qi in Eqsl~(26) and(2Y ) correspond to 



the label momenta of the collinear fields in Eqs. (16) and (17). The tree-level expressions for 



these hard functions for the various partonic channels are given in the Appendix. We note 
that the hard functions are known at the one-loop level 



In Eq. (23), large logarithms are summed via renormalization group (RG) equations for 
the hard, jet, and soft functions, while the PDFs are evolved via the standard DGLAP 



equations. In Eq. (23) large logarithms are summed by running hard function from the hard 



scale [in to the soft scale fis, while the jet function is run from the jet scale /ij to the soft scale 
fis- The soft function and the PDFs are evaluated at the soft scale /ig. The RG evolution 
equations for the hard and jet functions are given in the Appendix for completeness. 



5 Numerical Results 

5.1 Explicit form of the factorization formula at NLL 

In this section we give numerical results for the m e ff distribution for the pp — > 7 + 2 jets 
process. In particular, we present results at the NLL level of accuracy. We follow the 
conventions outlined in Table 1 of Ref. [25] for the counting of logarithms. At the NLL 



level of accuracy, the factorization formula in Eq. (23) can be computed using tree- level 
values for the hard, jet, and soft functions. The RG evolution is performed by running the 
hard function between the hard and soft scales fin and fi$- The jet function is run between 
the jet and soft scales /ij and fis- The soft function and the PDFs are evaluated at the 
soft scale. The RG evolution of the hard and jet functions at NLL accuracy requires the 
corresponding two-loop and one-loop cusp and non-cusp anomalous dimensions respectively. 
Finally, at NLL the LO PDFs are DGLAP-evolved using two-loop running of the strong 



coupling. After integrating Eq. (23) over f 2 in the range [0, f| ut ], the explicit form of the 
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Figure 3: The solid and dashed lines correspond to the m e ff distributions with NLL resum- 
mation and at LO, respectively. The qg initial state is shown in red while the q\q2 initial 
state, with $ representing any quark or anti-quark, is shown in blue. 



cross section at the NLL level of accuracy is given by 



dcr (pp -i> 7 + 2 jets) 
d$ qjr d$„ T . d$„ 



Ua 



dzS\ql+ql+q T )Tr \h^ s ^ h )S^ 



(2tt) 4 

xexp(-4(Cj 2 +Cj 4 )S(nj,fx s ) - Aj 2 (hj,h s ) - Aj 4 (/ij, /i s )) 



2/ij coshy 



Qj 4 



HhVs 



r(-u) 



fa{Xa)fb(Xb) 



(28) 



where the dependence of f| ut is implicit in the definition of z m i n . We have used the result of 



Eq. (71) to include the RG evolution of the jet functions in the above expression. We now 
explain the various parts of the above formula. The matrix H NLL (fis', ^h) denotes the tree- 
level hard function in a color- rotated basis with NLL RG evolution between /i# and fi$- in 
this color-rotated basis, the evolution of the hard function matrix elements is multiplicative, 
as seen in Eq. (43). The relation between the original hard-function matrix H obtained 



from matching QCD onto SCET operators is related to the function H in the rotated basis 
as shown in Eq. (44). Similarly, denotes the tree-level soft function in the same color- 
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rotated basis. The partonic momentum fractions x a ^ for the initial-state partons of flavor 
a, b are given by Eq. (25). The parameter u is defined as 



-2(Cj 2 +Cj 4 )A^j,^ s ), 



(29) 



and all the remaining quantities in Eq. (28 ) and on the right-hand side of Eq. (29 ) are defined 



in the Appendix. Note that for \in > [Ms, u < 0, so that the plus-prescription for the z — > 1 



limit has been dropped in Eq. (|28|). The rapidity y of the entire two-jet and photon system 

q% + q% + q z 



in Eq. (28) is given by 



y 



tanh 



q°j 2 + q°j 4 + q ' 



(30) 



is determined by 



max 



1 - -A(VA 2 + 4- A),rexp(2\y\) 



where 



A 



2 cosh y f 2 



cut 



TS 



(31) 



(32) 



For a fixed 



We choose f™ 1 to be of the order of the soft scale so that ff ut 
f| ut , in the machine-threshold limit r — > 1, z m % n in given by the second argument in Eq. (31 ) 
and is independent of f| ut . This is simply understood by the fact that in the threshold limit, 
the jets become extremely narrow and thus insensitive to f| ut . In other words, for a fixed 
f| ut the threshold condition eventually becomes a stronger constraint on the jet masses as 
compared to f| ut when r — > 1. For the numerical results, we choose f| ut to be large enough 
so that z m i n is independent of it. If we were to choose a small f| ut , logarithms of this quantity 
associated with vetoing additional jets would appear in the fixed-order result. Since T2 is 
related to 1 — z through the delta-function constraint in Eq. (23), our factorization formula 



also provides theoretical control over logarithms of this variable. This is seen explicitly in 
the NLL result above by the relation between f| ut and z min . We integrate over the phase 



space factors in Eq. (28) with the delta-function constraint S[m e ff — q j2 — qj 4 — q ] in order 



to generate numerical results for the m e ff distribution. 



5.2 Numerics for a 7 TeV LHC 

We begin by showing the m e ff distributions for the important partonic channels in Fig. [3j 
The dominant contributions come from the qiq2 and qg initial states, where the qi represent 
any quark or anti-quark. We note that the gg initial state contributes only at the percent 
level or less for the considered m e ff range. The dotted and solid curves correspond to the LO 
and NLL results respectively for the different partonic channels as described in the caption. 
Lower cuts of 100 GeV have been imposed on both jets and on the photon, and all final-state 
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Figure 4: The dependence of the ratio Rj on the choice of the jet scale [ij for various 
kinematic points. From the right to the left, are the curves used to determine the jet scales 
for m e ff = 1 TeV, 2 TeV and 4 TeV, respectively. 



particles are required to be separated by ARij > 0.4 with Ai2^ = (r]i — r/j) 2 + (0i — 02) 2 - 
CTEQ6 LO PDFs [31] have been used to produce these numerical results. In order to obtain 
the LO curves we choose the scales 



m = fij = us = y q T j 2 q T Jv (33) 

effectively removing the resummation. This dynamical scale is similar to scales found to 
reduce the effect of higher-order corrections in fixed-order calculations of W+multi-jet pro- 



duction [5J. For the NLL curves, we choose the hard scale as fin = Qj 2 Qj 4 - The soft 

scale is set by the see-saw relation /x 5 = y? H j ' The jet scale fij is determined numerically 
by minimizing the contribution of the logarithmic terms in the NLO jet function to the 
cross-section. In particular, we minimiz the quantity 

Rj = da NL °^-^(fi H = n s = /ij) 
da LO ({i H = hj = h s = \ q T h q T h ) 



with respect to fij. The numerator is defined as the contribution to the cross-section only 
from the logarithmic terms in the NLO jet functions. In Fig. [4], we show the dependence of 
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the ratio Rj on fij for several representative values of m e ff. We see that clearly identifiable 
minima occur for specific choices of \xj. As m e // increases the ratio fij/fiH becomes smaller, 
indicating a growing hierarchy between the hard and jet scales. As m e ff decreases, we see 
that the ratio fij/fiH becomes comparable to one, indicating that there is no longer a large 
hierarchy. In this case the effect of threshold resummation becomes smaller, and one can 
rely on fixed-order perturbation theory. For the numerical results presented in this section, 
the jet scale is determined numerically for each value of m e ff. We have tested that choosing 
other values of the hard scale, such as jiu = m eff/3, leads to negligible numerical differences. 



2.5 i— r 




Figure 5: A plot of the K-factor, defined as the ratio of the cross section with NLL resumma- 
tion over the LO result, for the important partonic channels. The curves from top to bottom 
correspond to the qg and qiq 2 initial states, where denotes any quark or anti-quark. 

From Fig. [3] we see that the effect of resummation becomes more important in the region 
of large m e ff, which is dominated by threshold kinematics. This is further illustrated in 
Fig. [5j where we show the K-factor for the same partonic channels as a function of m e ff. 
The K-factor is defined as the ratio of the cross-section with NLL resummation over the LO 
cross-section. Finally, in Fig. [6] we show the total K-factor and the uncertainties associated 
with scale variation. The solid red curve shows the total K-factor as a function of m e ff. The 
wider green band is obtained by setting //# = \lj — fig = \l and varying the scale \i in the 
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range {1/2, 2} around the central value \i = Jq T j <?J for the ratio 



Klo = ^ . (35 ) 



The narrower blue band is the result of scale variation of the NLL resumed cross-section. 
In this case, we vary the scales ^h,^j,Hs in the range {1/2,2} around their central values. 

The central value of the hard scale is taken as /i# = J q J 2 q J 4 . The determination of the jet 
and soft-scale central values was described earlier. More explicitly, the NLL scale variation 
band is determined by varying the scales Hh,^j,^s in the ratio 

da NLL {fi H ,fij,fi s ) x 
J<nll = , , (36) 

where again we have normalized with the NLL cross-section evaluated at the central scale 
choices fi c H , fij, fi c s . We see from Fig. [6]that scale variation uncertainty is significantly reduced 
when NLL resummation is included. 



6 Conclusions 

In this manuscript we have begun an investigation of the effect of threshold logarithms in 
gauge boson plus multi-jet production at the LHC. Such processes serve as backgrounds to 
supersymmetric particle production and to other forms of new physics that contain missing 
energy signatures. We have focused on pp — > 7 + 2 jet production, which is used as a 
calibration process for the production of missing energy in association with jets, as a first 
example. We have derived a factorization theorem using SCET that enables the resummation 
of large threshold logarithms. In our derivation we have used the global event shape N- 
jettiness to define the final-state jets. N-jettiness also allows for the resummation of large 
logarithms associated with vetoing additional jets, further extending the usefulness of our 
results. Numerical results comparing the leading-order m e ff distribution to the next-to- 
leading logarithmic resummed result have been presented. We find corrections ranging from 
50% to 100%. The logarithmic corrections increase with m e ff, potentially mimicking SUSY 
signatures, and should be accounted for in experimental analyses. 

Several future directions remain to be pursued. A next-to-leading order calculation of 
the soft function appearing in the factorization theorem would allow for an extension of 
the resummation accuracy to the NNLL level. The extent to which logarithmic corrections 
beyond NLO affect the high-m e // tail could be determined. The resummed result could then 
be combined with the fixed-order NLO calculation to better predict these backgrounds. An 
NLO calculation of the soft function would also allow the soft scale us to be determined by 
minimizing the contribution of the associated logarithmic terms to the cross section, as was 
done for the jet function. The degree to which the threshold region is dynamically enhanced 
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Figure 6: Plot showing the total K-factor defined as the ratio of the NLL and LO cross- 
section as a function of rn e ff- Also shown are the scale dependences of the LO and NLL 
results. The range of variation is defined in the text. 



in this process could then be determined. Our study defines jets via the N-jettiness event- 
shape variable, which is theoretically convenient since logarithms associated with jet vetoes 
can be controlled to all orders in perturbation theory. LHC experimental studies typically 
define jets via the anti-fey algorithm, and it is an open question as to the quantitative effect 
of this difference. Since N-jettiness jets are geometrically similar to anti-fey ones assuming 
the correct distance measure is chosen when defining jettiness [30], we expect the difference 
to be small. However, this point is worth further investigation. 

In the future we plan to extend these results to include W and Z production in asso- 
ciation with two and more jets, and combine the resummation with the known fixed-order 
results to provide a best prediction for use in experimental studies. We also plan to study 
the production of other color-neutral objects in association with jets, such as the Higgs bo- 
son. Since the Higgs is produced primarily through the gluon-gluon partonic channel, the 
threshold region can be enhanced even for a moderately energetic final state. In addition, 
experimental Higgs searches often divide the signal into exclusive jet bins, and control over 
jet- veto logarithms is needed. Our formalism handles both sources of logarithmic corrections, 
making for an interesting phenomenological application. 
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A RG Evolution 

Here we collect useful formulae for the evolution equations of the hard and jet functions. 
A.l Hard Function Evolution 



As seen in Eqs. (26) and (27), the hard functions H are determined via the Wilson coefficients 



C in Eqs. (18) and (19). Schematically, we can write the matching of the QCD operator O 



onto the SCET operators as 



O = 9id, 



(37) 



where the indices I run over the color structures shown in Eqs. (16) and (17). We have 



suppressed all spin structure. The hard function in Eqs. (26) and (27) can then be written 

as 



H 



u 



CjC 



(38) 



where we have again suppressed all spin-structure contractions on the right-hand side. 
The RG evolution equation in color space for the Wilson coefficients Cj is given by [26] 



dCV 
dlog/i 




(Ta ) 



u 



r cusp (a s ) log 



Ci = TijCj 



(39) 



which is valid at least up to two loops for N external massless colored particles. Here, 
s a b = 2<r a feg a • qb + iO with a a b = 1 if both q a and qb are incoming or outgoing, and a a b = — 1 
otherwise. q a ,b denote the label momenta on the external fields on the SCET operators (see 



Eqs. (16) and (17)), and the indices a or b run over all the external fields. The logarithm of 



s a b can then be written as 

l0g(-S a6 ) = log \s ab \- A af) 27T , (40) 

with A a b = 1 for a and b both incoming or outgoing and otherwise. The action of the color 



matrices T of Eq. ( 39 ) on the collinear quark and gluon fields is defined as [27] 



= -t A £ b 5, 



abi 



&abibk A , 
5abAfrf CAB . 



(41) 
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We note that in the evolution of the SCET operators, the Tjj on the right-hand side of 
Eq. (39) acts on the 9 from the right: OiYij. 

For the ab — > cd + 7 processes considered here, with a, b, c, d denoting colored patrons, 
the anomalous dimension Fjj defined in Eq. (39) can be explicitly written as 

2 4 

(c 3 + c 4 ) 5/jr cusp log h V] i a bu 

— S34 ~~ 
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(ci + c 2 ) 5/jr cusp log 

z — S12 



1 



rpArpA 
J l J 3 1 



cusp 



log — — + log 



2 

-S34 



-S13 -Sl3 



+ 



rpArpA 
1 1 1 A 1 



cusp 



log 



-512 
-Sl4 



log 



-S34 



cusp 



, — s 12 , — S34 

log h log 

-S 2 3 -523 



rpArpA 
1 2 1 i I 



cusp 



log — ^ + log 

-S24 



Su 

— S34 

-S24 



/J 
(42) 

where we have used the fact that — due to color-charge conservation and the relation 

T a ■ T a = c a , where c a = Cp for quarks and c a = Ca for gluons. The RG evolution equation 
of the hard coefficient Hjj can then be determined by Eqs. (38) and (39). The result is given 
by 

Hij(ij) = Hu(hh) exp(2c H S(nH,p) - 2A H (p H , (jl)) 



xexp -Ay(p H , fi) —log 



-si 2 )(-s 34 ) 



f4i 



+ 2c r log 


— S24 




-S13 



+ a, + a: 



(43) 

where H is related to H by the transformation 

H = PHP ] . (44) 

The transformation by the matrix P corresponds to a rotation in color space that diagonalizes 
the matrix V on the right-hand side of Eq. (39). In this basis, the evolution of the matrix 
elements Hjj in color space is multiplicative, as seen in Eq. (43). The A/ in Eq. (43) are 
just the eigenvalues of the matrix T, and c r is a constant. The Aj and c r constants depend 



on the partonic process. Finally, the remaining quantities in Eq. (43) are defined as 

' QbM r cus P (a) d «' 



da 



(3(a) 

Fpusrifa) 



i cusp 1 

° M dQ E^ 



a s (pii) 



0(a>) 



(45) 



To NLL accuracy, these quantities take the form 



Air 



A0l \a s (/ii 
To 



2/?o 



logr + 



1 - - 

r 

Air 



logr + 



£1 

r 



fa 

0o 



ri 


01 


r 


00 


(r - 





[1 —r + logr) + 



I 
2/3, 



r log 2 r ] 



(46) 
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where 



7° 



a s (/i f )/a s (/ii 



4, r 1 = 4 



67 TX< 

"9 3 



C A 



34 



ll r 2 r 



9 
2 



2 : 



7<? 

7 9 



-3CV — for quarks, 

4:71 



a*. 



for gluons. 

4tT 



(47) 



A. 1.1 Four External Quarks 



Processes with four external quarks in Eq.(14) are mediated by the operators 



Ox 

0-2 



616 616, 



(48) 



where we have suppressed the Dirac structure. With this choice of labels for the collinear 
fields, which is same as in Eq. (16)),we have T\ ■ T 3 — T 2 • T 4 , T x ■ T 4 = T 2 ■ T 3 where 
T a -T b = J2 A TfT b A and 



{Ti-T 3 ) xl 



1 

"C~A 

C A 



(Ti ■ T, 



3j21 



Cf 
2C> 



(7\ ■ T, 



3J12 



(7\ • Tj 



3; 22 



2 + C A 



(Ti • T, 



4j21 



2C7/ 



4H2 



0. 

{Ti ■ T 4 ) 22 = . 

(49) 



Using these results in Eq. (42), the anomalous dimension matrix takes the form 

-S12X-S34) 



1 j 



+ 



(4CV)r cusp - log 

T x ■ T 3 



T 

2 cusp 



2 log 



y~ S 12)( _ S34) \ Tl • T 4 



(50) 



+ 



r 

2 x cusp 



" s 13)( — s 24 

This corresponds to the values 

c H = 4C F , ^ 7 a = 4 7g , c r = 



2 log 



-Sl2)(-S34) 
-S23)(-Sl4) 



(51) 



in Eq. ( 43 ) . The eigenvalues A/ of the matrix T are given by 

\± = ^{U-T)-±U±yjuT+ \C\{T - Uf 



(52) 



19 



where we have defined 

TT 1, (-auK-ggg) rp 1, (-gl2)(-g 34 ) , -x 

2 (-S13X-S24) 2 (-Si3)(-s 24 ) 

The transformation matrix P which diagonalizes V is given by 

p= \ A + p\ p-^^^(% v -%sp 

"' 1 A_ gi7 y ' AXCpU V -A- A+ 

A. 1.2 Two external quarks and gluons 

The processes with two external quarks and gluons are mediated by the operators 

0! = &t a H as ^A ai A a '\ 

2 = &t a H ai £ 4 A ai A a3 , 

3 = {= 2 5 aia3 £ 4 A ai A a3 , (54) 

where we have again suppressed the Dirac structure. With this choice of labels for the 
collinear fields we have T\ ■ T 4 = T 2 ■ T 3 and 

(2~i ■ 23)11 = (2\ • T 3 ) 22 = — — Cp , (2\ • 23)31 = (Ti • T 3 ) 32 = — - , (2i • T 3 ) 33 = — Ca 

zl-a 4 

(Ti-T 4 ) 3 i = j, (Ti-T 4 ) 22 = -^, (T 1 -T 4 ) 13 = -(T 1 ^4)23 = 1, 

(T 2 -T 4 ) n = (T 2 -T 4 ) 22 = -^, (T 2 -T 4 ) 3 i = (T 2 -T 4 ) 32 = -j, (T 2 ■ T 4 ) 33 = -C F . (55) 

ZCa 4 

We have only listed the non-vanishing elements. The anomalous dimension matrix is given 
by 



u = \ ^(2C A + 2C F )T cusp - log v 1Z ; V 4 ^ + 2 lq + 2 lg ) 5 U (56) 



. C F — Ca\ r , -s 24 

+ ( — z — r cusp log — d/j 

^ / — Sl3 

2\ • T 3 T 2 -T± \ (-gigK-ggp , 2 1 i-T 4 / (-si 2 )(-s 34 ) 
+ 1 ( — o 1 n — ) r cus P log t TT r H — r cusp ( 2 log 



2 2 y — (-Si3)(-a M ) 2 — V (-^i4)(-s 23 

where we have used (2\ ■ T 3 — T 2 ■ T 4 ) 7J = (Cp — Ca)5ij- This corresponds to the values 

2C A + 2CV, ^7 a = 2 75 + 2 7g , c r = ^(C F - C A ) , (57) 



7 j 



ch 



in Eq. (43). The eigenvalues of V and the rotation matrix P can be obtained from |26j. To 



maintain some semblance of brevity we do not provide explicit expressions for them here. 
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B Tree-level expressions for the hard functions 



The hard functions Hjj that appear in the factorization formula of Eq. (23) take the forms 
shown in Eqs. (26) and (27). In this section, we collect the tree- level expressions for these 
hard functions for the different partonic channels that contribute to pp — > 7 + 2 jets. The 
partonic channels in Eqs. (14) and (15) are mediated by the SCET operators in Eqs. (16) 
and (17), respectively. The collinear fields in these operators are labeled by indices {1,2,3,4}. 
In our convention we fix the basis of operators with labels as in Eqs. (48) and ( |54| ). The 
momenta of the particles in the partonic process are chosen to correspond to a permutation 
of the label momenta {gi, q 2 , (ft, #4} of the fields. Thus, depending on the partonic process, 
these label momenta will correspond to a particular combination of incoming and outgoing 
momenta. The hard coefficients are then given in terms of these label momenta that are 
assigned in the partonic process. This allows for a consistency in the momentum assignments 
in the matrix element calculations and those that appear in the RG evolution equations. By 
convention, we will denote the outgoing photon momentum by q$ in the following. We 
note that all hard function have been checked to numerically agree with those found using 
Madgraph 



B.l q{qi)q{qA) ->> QMQfehfe) 

Here we give the hard matching coefficient for the process q(qi)q(qi) — > Q{.q2)Q{.qz)l{.l5) 
with different quark flavors q 7^ Q. The matching coefficient is given by 



C F 



2 



Cf <±_ 
C\ C\ 



where the four-vector I is given by 



7 = 2 -e q h e q e Q h e Q — . (5£ 

\ S45 S15 S35 S25 / 

e q ,Q denote the electric charges of the quarks appearing in the scattering process. 
B.2 q(qi)Q(q A ) q{q2)Q{qi)l{q^) 

In this channel with q 7^ Q, the hard function matrix is given by 



4ATf ' v s > V s 12S34 J \ 

where the four-vector I is given by 

94 q\ qs . q2 \ , KO , 

1 — + e g e Q — + e q — J. (59) 

S45 S15 S35 s 25 / 
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B.3 q(qi)q(qi) -> 9(92)9(93)7(95) 

In this channel with identical quark flavors, the hard function is given by 

H = H a + H b + H c , 



where 



H* = ^ (4^ (-4^) ( <±±±£±£ ) P ( J I ) 



The four-vector / is given by 



S14S23S12S34 



1 

cl 

C F 

cl 


C F ' 


) /21 


V 2C A 



. , . , 94 ,91 93 . 92 

I = 2 I -e g h e g e g h e g — 

S45 Sl5 s 35 S25 



B.4 q{qi)Q{q 3 ) -> 9(92)^(94)7(95) 

In this channel for 9 7^ Q the hard function is given by 



H = — (4tt) 3 f-4aa 2 ) ( s i3 + s 24 + gi4 + 4s \ ( 1 



where the four-vector / is given by 



, , 9i , 9s . 92 . 94 
/ = 2 I h e Q h e 9 h e Q — 

Sl5 s 35 $25 S45 



B.5 9(91)9(93) 9(92)9(94)7(95) 

In this channel with identical quarks the hard function is given by 

H = H a + H b + H c , 



where 



and 



H a = — (4tt) 3 (-Aaa 2 ) ( S ' 3 + ^ 4 + S ' 4 + ^ 3 ^ I 2 ( 1 ° 
4iV c 21 J ^ s; V S12S34 J V 



1 C F 



ff b = 7^(47r) 3 (-4aa a 2 ) f S ' 2 + §34 + S ' 3 + 4 A I 2 [ ^ g< 
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AN? 



(4tt) 3 (-4cm s 2 ) (s 2 13 + s 



24) 



S13S24 — S12S34 — S14S23 
S14S23S12S34 



_ J_ Il£L 
Ca 2Ca 

2C A U 



The four-vector / is given by 

1 = 2 



? Sl 5 



Q3_ 
'S35 



's 2 5 



(U_ 

'S45 



B.6 g(q u ai)g{q 3 , a 3 ) ->• q{q2)q{qA)l{q^) 

For this gluon-initiated process, we find the hard function 

1 



H 



(I A \3 2\ 1 a 2 ^i=l,3,5 S *2Si4 (sf 2 + ^4) 

« 4 *> 16e * au^; 



x 



S14S23 |(si3«24 - S14S23 - S12S34) 

|(Sl3S24 - S14S23 - S12S34) S12S34 















B.7 q{qA)q{q2) g(qi, a{)g(q 3 , a 3 )-f(q 5 ) 
For this scattering process we find the hard function 



H 



1 (( 47 r) 3 ^)l6eJ E ™ Sl2St4(4 + 4) 



4iV c 2 



S 13 lL=l,3,5 S *2Si4 



S14S23 |(si3«24 - S14S23 - S12S34) 

> I 1(513^24 - S14S23 - S12S34) S12S34 





(64) 



B.8 q(q±)g(qi) ->• gfe)0fehte) 

For the gg channel we obtain the hard function 

IT 1 ( (a \3 2\ -.a 2 ^i=l,3,5 S *2Sj 4 (sf 2 + S i4) 



X 



S 14S23 ^(Sl3 s 24 — s 14 s 23 — s 12 s 34) 

2 ( s 13 s 24 — S 14S23 — s 12 s 34) Sl2 s 34 





It takes exactly the same form for other qg channels, such as q(2)g{l) — > 9(4)^(3)7(5). 
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C SCET jet and soft functions 



In this section, we present the definitions of both the jet functions and the soft function used 
in our approach. The jet functions defined in the framework of SCET have been known for 
some time up to the two-loop order [TBI H22 EH ES] • We refer the readers to the operator 
definition of the jet functions therein and we note that the jet functions J(s) used in our 
current work are normalized to S(s) at LO. The soft function showed here characterizing the 
restricted soft radiation inside and outside the jet region near threshold is new. 
The operator defintion of the soft function is 



Sji (k° out , t J n ] s ) = (0| O s J 6 (k° out -k°)l[6 (rj*, - fjf ) 0\ |0> , 



(65) 



where k° and t$ are the operators which act on the final ultrasoft states to project out the 
related observables. O s is a collective of ultrasoft Wilson lines Y and color structures T 
appearing in the SCET operators in Eqs. (16) and (17), which formally can be written as 

O" T Y^TY.YlTY, . (66) 



The subscripts of Y's correspond to the subscripts of the fields in each SCET operator. 

In this manuscript, the jet functions are run from the jet scale to the soft scale via the 
RG evolution equation 



d log n 



ds' 



or v — 

^^a 1 cusp 2 



+ 7 a As - s>) 



J a {s',fi) . 



(67) 



The index a runs over {q, g} corresponding to quark and gluon jet functions respectively, 
C q = Cf and C„ = Ca- To leading order, 



a-.s 



2/3, 



! 47r' 



(68) 



In position space with F(y) = f dsexp(—isy)F(s), this equation becomes 

^y±= (2C a r cusp iog( w V-) + 7 }) J a (v). 

(1 lOg /! 

The solution is 

J a (y,fi) = exp(-AC a S(fij,fi) - Afoj,?)) (iy^y 2Ca ^ M 



Ja(y,lij) 



(69) 



(70) 



Here Aj is obtained by replacing r cusp by 7} in Ay defined in Eq. (46). For the product of 
jet functions that appears in the factorization theorem, we have 



J a J b = exp(-4(C a + C b )S -A a j- A b j) (ly^e**) 



-2(C a +C 6 )A 7 



Ja{y,nj)Jb{y^j) 



(71) 
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The NLO jet functions are used to determine the jet scale in our current work and are 
given by 

Us) = 6(s) + ^ ((7 - * 2 Ms) - jM s lf) + ^ s ^ 2 )) > 

Us) = 5(a) + g ((Q - vr 2 ) C A + 6(a) - ^£o(s/f) + ^A(.//i 2 ) ) • 

(72) 

For NLL accuracy we need only the LO soft function which takes the form 

S?l = SjiS{k° out ) . (73) 

Ji 

where for the four-quark channels the color matrix S is diagonal and reads 




while for the two-quarks and two-gluon configurations the color matrix S is 

^A^F ~2 L ' F °A^F \ 
— \Cp CaC\ CaCf J • 

CaCf CaCf ^CpC\ J 
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